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1.0  Optoacoustic  Studies  of  Liquids 


1.1  Introduction 

Investigations  of  optoacoustic  generation  in  liquids  are  well 
12  3 

documented.  ’  ’  Most  of  these  studies  were  carried  out  using  low 
optical  energies  and  small  optical  absorption  coefficients.  A  small 
optical  absorption  coefficient  makes  the  symmetry  of  the  problem 
simple;  the  mathematics  tractable.  In  the  low  optical  energy  regime, 
the  acoustic  generation  mechanism  is  thermal  expansion.  This  is  the 
basis  for  linear  optoacoustic  theory.  Nonlinear  effects  occur  when 
the  liquid  experiences  multiphoton  absorption  or  a  phase  transition. 

On  a  molecular  level,  the  optoacoustic  effect  can  be  broken 
down  into  several  different  steps.  The  process  of  converting  optical 
energy  to  acoustical  energy  begins  with  the  absorption  of  photons  by 
the  molecules  of  the  fluid.  This  absorption  leaves  the  molecules  in  an 
excited  state.  In  the  present  case,  a  nitrogen  laser,  used  as  the 
energy  source,  produces  electronic  excitation  of  the  molecules. 
Dissipation  of  this  energy  can  occur  through  several  mechanisms 
including;  photochemical  processes,  radiative  decay,  and  nonradiative 
decay.  In  nonradiative  decay,  the  electronic  excitation  energy  is 
transferred  to  translational,  rotational  or  vibrational  states  of  the 
molecule.  Through  interactions  with  other  molecules  this  rotational 
and  vibrational  energy  can  be  converted  into  translational  energy 
creating  a  temperature  rise  in  the  fluid.  In  our  experiment,  following 
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each  pulse  from  the  laser,  the  fluid  is  allowed  to  cool  to  an  ambient 
temperature.  This  alternate  heating  and  cooling  creates  a  pressure 
variation  which  is  observed  as  acoustic  pulses. 

4  5  6 

Present  optoacoustic  theories  ’  ’  assume  that  the  energy 
transfer  processes  discussed  above  all  take  place  over  a  period  of 
time  which  is  very  short  compared  to  other  factors  that  effect  the 
observed  pressure  profile;  such  as,  the  length  of  the  laser  pulse  and 
its  diameter  as  it  enters  the  fluid. 

The  process  of  energy  transfer  and  the  rate  at  which  it  occurs 
varies  from  one  liquid  to  another  because  the  energy  levels  and  rate 
of  energy  transfer  depend  on  the  structure  of  the  molecules.  If  the 
rate  of  production  of  translational  energy  is  slow  due  to  internal 
relaxation  processes,  the  shape  of  the  acoustic  pulse  should  be 
governed  by  these  relaxation  processes.  The  experimental  system 
must  be  capable  of  resolving  very  short  acoustic  pulses  in  order  to 
resolve  events  occurring  on  this  time  scale.  Otherwise,  the  pressure 
profile  is  dominated  by  geometrical  factors  of  the  experiment. 

Carbon  disulfide  (CS2)  has  been  chosen  as  one  subject  of  study 
because  previous  ultrasonic  absorption  measurements7  have  already 
shown  that  the  transfer  of  energy  between  translational  and 
vibrational  states  is  relatively  slow  (compared,  for  example,  to  H20). 

Mr.  Charles  Thompson  completed  his  M.S.  work  in  December  of  1987 
which  was  a  study  of  energy  transfer  in  CS2.  His  Thesis  was  issued 
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as  a  Technical  Report. 

1.2  Experimental  Configuration 

The  experimential  setup  will  be  described  here  for 
completeness  even  through  a  similar  description  has  been  provided 
in  previous  reports. 

A  block  description  of  the  experimental  technique  is  shown  in  Figure 
1.1.  This  technique  utilizes  two  lasers.  A  1.25  MW  nitrogen  laser 
emitting  800  ps  pulses  in  the  UV  (337nm)  is  used  as  the  optoacoustic 
source.  A  3mW  continuous  He-Ne  laser  runs  parallel  to  the  excitation 
zone.  It  is  then  focused  on  a  photodiode  which  is  coupled  to  a 
digitizing  oscilloscope.  The  probe  laser  responds  to  refractive  index 
gradients  which  are  produced  as  a  result  of  the  acoustic  wave.  These 
gradients  cause  small  deflections  of  the  probe  beam  which  are 
detected  by  the  photodiode.  The  observed  signal  is  proportional  to 
the  pressure  gradient.  Comparison  to  theory  can  be  achieved  by 
either  integrating  the  observed  signal  or  by  differentiating  the 
predicted  pressure.  Both  techniques  have  been  used;  however,  a 
comparison  to  the  absolute  magnitude  of  the  acoustic  pressure  can  be 
achieved  as  well  if  the  theory  is  compared  to  the  output  of  a  PVDF 
hydrophone  as  shown  in  Figure  1.2.  Given  the  inaccuracies  in  beam 
size  estimates,  the  agreement  between  theory  and  experiment  is 
considered  acceptable. 
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Energy  delivered  to  the  fluid  is  varied  by  placing  glass  slides  in 
the  beam  path  of  the  UV  laser.  N  number  of  glass  slides  (-.23  mm 
thick)  are  supported  by  two  small  pieces  of  plexiglas.  The  centers  of 
the  plexiglas  have  been  removed.  To  provide  for  accurate  energy 
measurements,  twenty  sets  of  attenuators  corresponding  to  a 
different  number  of  glass  slides  were  made.  The  energy 

transmission  for  each  attenuation  was  measured  with  an  energy 
meter.  The  precision  of  the  energy  meter  is  ±.  lpJ.  Figure  1.3  shows 
measured  energy  delivered  to  the  test  cell  as  a  function  of  the 
number  of  glass  slides. 

The  energy  density  in  the  fluid  can  also  be  increased  by 
increasing  the  optical  absorption  coefficient  a.  This  is  achieved  by 
either  adding  dye  to  the  test  liquid  or  by  choosing  a  liquid  that 
strongly,  absorbs  UV  radiation  (CS2). 

Determination  of  the  Excitation  Beam  Size 

The  dimensions  of  the  excitation  beam  are  among  the  most 
important  parameters  involved  in  theoretical  predictions.  An  accurate 
measurement  of  the  spatial  beam  profile  is  vital  for  comparison  to 
theory. 

The  Nitrogen  laser  used  emits  an  800  ps  UV  pulse,  rectangular  in 
shape.  The  ratio  of  the  two  sides  of  the  rectangular  beam  can  be 
determined  by  passing  a  small  piece  of  photographic  film  through  the 


4 


beam  in  the  vicinity  of  the  focal  point.  The  size  of  the  exposed  portion 
of  a  film  measured  under  a  microscope  was  found  to  be  150  pm  x  550 
pm.  Although  the  absolute  dimensions  will  depend  upon  the  shape  of 
the  beam  wings  and  exposure  time,  the  width  to  height  ratio  550/150 
should  remain  constant. 

The  time  duration  of  the  acoustic  pulse  is  primarily  determined 
by  the  width  of  the  short  dimension  of  the  excitation  pulse.  By 
comparing  an  experimental  waveform  to  theory  one  can  estimate  the 
small  dimension  of  the  rectangular  pulse.  For  this  measurement,  it  is 
important  to  use  a  liquid  with  a  fast  relaxation  time  to  insure  that  the 
shape  of  the  acoustic  pulse  is  not  effected  by  relaxation  effects. 
Fluorescence  may  also  cause  problems  in  beam  size  measurements  due 
to  radiative  transport  outside  the  initial  excitation  zone. 

Alcohols  and  water  are  fast  relaxing  fluids  but  are  not  good 
absorbers  of  UV  radiation  (337  nm).  One  way  to  increase  the 
absorption  rate  is  to  add  dye  to  the  fluid,  however  the  dye  presents  a 
problem  with  fluorescence.  The  liquid  chosen  for  the  measurements 
presented  here  is  a  clean  brand  of  Scotch.  It  is  alcohol  based,  a 
reasonably  good  absorber,  relaxes  rapidly,  and  fluorescences  only 
weakly.  This  fluid  provided  the  shortest  optoacoustic  signal.  The 
comparison  of  experimental  and  theoretical  waveforms  are  shown  in 
Figure  1.4.  The  radius  required  to  match  the  waveforms  was  45  pm. 
This  suggests  a  rectangular  profile  of  90  pm  x  330  pm.  Using  this 
technique  the  values  given  should  be  accurate  to  within  £  10  pm . 
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1.3  Microscopic  Interpretation  of  Optoacoustics  in  Fluids 


At  the  UV  wavelengths  employed  in  this  experiment,  the  laser 
photons  whicK  are  absorbed  result  in  electronic  excitation  of  the 
recipient  molecules.  This  electronic  energy  is  then  rapidly  converted 
to  vibrational  a^d  translational  energy  which  is  manifested 
macroscopically  as  thermal  expansion.  We  conclude  that  the 
electronic  excitation  does  not  lead  directly  to  thermal  expansion  by 
comparison  of  the  UV  results  to  measurements  using  IR  lasers.  The 
IR  lasers  produce  photons  with  insufficient  energy  to  cause  electronic 
excitation  yet  the  efficiency  at  which  the  electromagnetic  energy  is 
converted  to  thermal  expansion  is  identical  to  that  for  the  UV 
measurements  (at  least  in  the  linear  regime).  Since  the  observed 
thermal  expansion  depends  only  upon  the  total  energy  deposited  in 
the  fluid  and  not  on  the  optical  wavelength,  we  conclude  that  it  is 
vibrational  and/or  translational  modes  which  give  rise  to  expansion. 

Further  insight  into  the  critical  energy  flow  path  can  be  gained 
from  Thompson's8  results.  Specifically,  in  the  case  of  CS2,  he  found 
that  following  deposition  of  energy  by  the  laser,  there  was  a  time  lag 
for  resultant  thermal  expansion  of  about  90  nsec  which  seemed  to  be 
associated  with  transfer  of  energy  from  an  electronic  (initial)  state  to 
translation  (or  vibration)  as  well  as  a  time  constant  of  1.5  nsec  for 
vibrational  relaxation.  These  results  indicate  that  when  energy  is 
added  to  the  fluid,  that  energy  must  undergo  internal  conversion 
processes  prior  to  causing  thermal  expansion.  Further,  since  the 
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laser  used  has  a  pulse  length  of  0.8  nsec,  essentially  all  the  laser 
energy  deposited  resides  in  internal  states  at  the  end  of  the  pulse. 

The  observation  that  deposited  energy  resides  in  excited  states 
for  relatively  long  time  periods  gives  rise  to  some  interesting 
consequences.  Specifically,  for  a  large  energy  flux,  this  raises  the 
possibility  of  depositing  photons  in  a  grouping  of  molecules  more 
rapidly  than  the  energy  is  redistributed  in  the  fluid. 

Refer  to  Figure  1.5.  Up  to  an  input  energy  of  about  80  pj,  the 
optoacoustic  signal  strength  varies  linearly  with  input  energy.  The 
slope  of  the  curve  is  consistent  with  linear  thermal  expansion.  At  80 
pJ,  there  appears  to  be  a  plateau  followed  by  a  range  of  energies 
where  the  signal  strength  increases  more  rapidly.  At  about  160  p.J, 
there  is  another  plateau  followed  by  an  even  steeper  rise  in 
optoacoustic  signal  strength  with  input  energy.  When  the  input 
energy  reaches  210  pJ,  the  slope  begins  to  decrease.  There  is  a 
further  decrease  at  320  (iJ  returning  the  slope  to  near  its  low  signal 
value.  It  is  interesting  to  note  that  at  each  of  the  two  plateaus,  the 
optoacoustic  conversion  efficiency  (the  slope  of  the  curve)  increases 
by  approximately  50%. 

Higashigaki,  et  al.9  has  carried  out  investigations  of  molecular 
reorientations  of  carbon  disulphied  molecules  using  a  light  scattering 
technique.  The  reorientation  of  the  CS2  molecular  axis  involves 
displacements  of  other  molecules.  He  concludes  that  the  activation 
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energy  required  to  reorient  a  CS2  molecule  is  1.58  kcal/mol. 
Considering  the  size  of  our  excitation  volume,  we  find  that  there  are 
1.33  x  10' 8  moles  of  CS2.  This  gives  an  activation  energy  of  2.1  x 
10'5  cal  or  88  pJ.  Note  that  this  occurs  very  near  the  first  plateau 
suggesting  that  these  reorientations  may  contribute  to  the  increase  in 
optoacoustic  conversion  efficiency. 

The  theory  from  which  the  beam  size  parameters  were 
deduced10  assumes  a  Gaussian  spatial  profile.  Although  a  radial 
profile  is  assumed  there,  for  purposes  of  energy  disposition 
calculation,  we  will  adopt  a  rectangular  Gaussian  described  by 


I  =  Io  e-*2/x02  e-y2/y0  2  e'“ 


1 

I 

I 

I 

I 


where  x0  is  estimated  to  be  90  |X,  y0  is  estimated  to  be  330  p.,  and 
1/a  has  been  measured  to  be  27  p..  From  the  temporal  response  of 
the  optoacoustic  waveform,  we  conclude  that  the  energy  contribution 
is  substantially  constrained  to  a  volume  x0  x  y0  x  1/a.  Some  of  this 
region  will  absorb  more  energy,  other  parts  less.  We  will  comment 
on  this  fact  later. 

Based  upon  the  above,  as  a  starting  point,  we  will  assume  that 
the  beam  has  an  effective  cross  section  of  90  p.  x  330  p  giving  an 
area  of  3  x  lO'^m^.  We  will  further  assume  that  the  pulse  length  is 
800  psec  (manufacturer's  specification)  so  the  input  intensity  at  the 
first  plateau  (80  pJ)  is  3.3  x  lO^w/m^  .  In  accordance  with  Beer’s 
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Law, 


I  =  I0e-«z 

where  I0  is  the  input  intensity,  a  is  the  absorption  coefficient,  and  z 
is  the  depth  into  the  fluid.  Differentiating, 

dl/dz=  -aI0e-(Xz 

gives  the  energy  absorbed/unit  volume-unit  time,  a  distance  z  into 
the  test  cell.  At  the  point  where  the  excitation  beam  intensity  is 
down  to  1/e  of  the  initial  value,  az  =  1  and  dl/dz  =  -al0/e.  Taking  a 
=  3.7  x  lO^m'l  (penetration  depth  of  27  p.m),  dl/dz  =  4.5  x 
lO^j/mS-s  at  z  =  1/a  and  1.2  x  10*7  J/m3  -s  at  z  =  0. 

The  results  of  this  energy  deposition  can  be  viewed 
macroscopically  and  microscopically.  Microscopically,  the  energy 
absorbed  can  be  viewed  in  terms  of  the  number  of  photons.  Taking  X 
=  3.4  x  10'7m,  E  =  hc/X  gives  an  energy  of  5.9  x  10_19J/photon  so  the 
fluid  is  absorbing  between  .76  and  2.1  x  1035photons/s-m3.  For 

comparison,  at  20°C,  the  density  of  CS2  is  about  1.3  g/cm3  and 
since  each  CS2  molecule  has  a  mass  of  1.3  x  10_22g,  there  are 
approximately  1022  molecules/cm3  or  1.03  x  1028  molecules/m3.  If 
the  excitation  laser  operated  in  the  steady  state,  each  molecule  would 
absorb  1-2  x  107  photons  but  since  the  pulse  is  only  800  psec  in 
duration,  each  molecule  absorbs  only  .02  photons  during  the  laser 
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pulse.  Another  way  to  view  this  is  that  if  the  photon  density  exceeds 
about  0.02  photons/molecule,  a  plateau  in  optoacoustic  conversion 
energy  is  reached  —  something  different  begins  to  happen. 

Next,  we  will  adopt  a  macroscopic  view.  In  the  defined 
excitation  volume,  the  energy  absorbed  near  the  first  plateau  is 

about  1.0  x  108j/m3.  The  specific  heat  of  CS2  liquid  near  room 

temperature  is  76  J/g-mol-K12.  The  density  is  1.3  g/cm3  or  0.017  g- 
mol/cm^  so  Cp(/)  is  1.3  x  lO^J/m^-K  resulting  in  a  temperature  rise 
(ignoring  temperature  dependence  of  density  and  Cp(/))  of  77*  K . 
The  boiling  point  of  CS2  is  46* C.  Assuming  a  starting  temperature 
of  23*C,  the  plateau  at  80  pJ  is  in  excess  of  the  energy  necessary  to 
raise  the  fluid  to  its  boiling  point  throughout  the  focal  volume. 

Two  minor  details  should  be  noted  prior  to  continuing  this 

argument.  First,  the  excitation  laser  intensity  is  greatest  next  to  the 
quartz  surface  at  the  beam  entry  location.  The  fluid  is  in  contact 
with  a  surface  which  should  prohibit  supersaturation.  Second, 

microbubbles  are  not  observed  until  the  input  energy  reaches  about 
120  pJ.  At  an  energy  required  to  produce  bubbles,  no  unusual 
dependence  of  optoacoustic  signal  strength  with  input  energy  is 
observed. 

1 3 

The  heat  of  vaporization  of  CS2  is  89.56  g  cal/g  or  4.9  x 
108j/m3.  The  first  2.7  x  lO^J/m^  is  required  to  reach  the  boiling 
point  (prior  to  the  80  pJ  plateau).  The  next  4.9  x  lO^J/m  3  is 


required  to  provide  the  heat  of  vaporization  at  z  =  0,  the  total  energy 
required  to  vaporize  the  fluid  would  be  400  pJ.  Note,  however,  that 
between  the  80  p.  J  and  160pJ  plateaus  small  bubbles  begin  to 
appear  and  the  slope  of  the  curve  increases.  This  suggests  that  in  the 
focal  volume  where  the  laser  intensity  is  at  its  highest,  local 
vaporization  occurs  at  lower  energies 

From  a  macroscopic  point  of  view,  the  signal  strength  versus 
input  energy  is  now  clear  at  least  up  to  the  point  near  210  pJ  where 
the  slope  beings  to  decrease.  Specifically,  at  low  input  energies,  the 
optoacoustic  amplitude  increases  linearly  with  input  energy  up  to  80 
pj.  At  this  point,  sufficient  energy  is  deposited  to  bring  the  fluid  to 
the  boiling  point.  Further  increases  in  input  energy  serve  only  to 
provide  heat  of  vaporization.  Once  some  of  the  fluid  is  vaporized,  the 
coefficient  of  thermal  expansion  becomes  dominated  by  those  gas 
phase  molecules  present.  The  gas  phase  has  an  expansion  coefficient 
(assuming  an  ideal  gas)  of  3.1  x  10'3/K  while  the  liquid  phase  has  an 
expansion  coefficient  of  1.2  x  10'^/K.  Referring  to  Figure  1.5,  the 
slope  of  the  curve  at  low  energy  is,  indeed,  0.4  times  that  above  180 
pj.  Between  these  two  slopes,  there  is  an  intermediate  region  where 
part  of  the  energy  is  being  used  to  do  expansion  work  but  most  is 
going  into  the  heat  of  vaporization. 

Now  that  the  general  features  of  the  experimental  observation 
have  been  explained  from  a  macroscopic  viewpoint,  again  consider 
what  is  happening  at  the  molecular  level.  To  begin  this  description, 
one  more  important  experimental  observation  needs  to  be  noted. 


A 
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Specifically,  as  the  energy  deposited  in  the  fluid  increases,  the 
optoacoustic  waveform  does  not  change  in  shape.  Since  the  shape  is 
controlled  by  internal  relaxation  processes,  this  observation  proves 
that  energy  must  appear  in  translation  and  vibration  prior  to  an 
observable  expansion  or  change  in  state.  We  can  and  will  assume, 
therefore,  that  all  processes  at  the  molecular  level  start  after  energy 
is  transferred  to  v-t  states. 

We  will  view  the  liquid  as  a  dense  gas  or  imperfect  lattice 
characterized  by  some  radial  distribution  function  g(R)  which  gives 
the  probability  of  finding  a  second  molecule  a  distance  R  from  a 
reference  molecule.  For  a  perfect  lattice,  g(R)  would  be  periodic 
along  any  axis;  for  a  gas  at  low  density,  it  would  be  zero  out  to  the 
molecular  diameter  and  uniform  beyond.  Since  we  envision  the 
liquid  as  an  intermediate  state,  g(R)  should  be  a  damped  periodic 
function  which  decays  toward  a  constant  value  more  rapidly  with 
increasing  temperature  until,  at  the  boiling  point,  g(R)  approaches 
the  gas  phase  value.  Indeed,  microscopic  description  of  the  liquid 
state  typically  involves  expressions  for  g(R)  in  terms  of 
intermolecular  potentials. 

Equation  of  State 

We  have  shown  that  for  large  deposition  energies,  the  liquid  in 
the  excitation  volume  may  undergo  a  change  in  phase.  Heating  of  the 
volume  will  change  the  thermodynamic  properties  of  the  fluid  which 
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enter  into  subsequent  calculation. 

An  equation  relating  the  volume,  temperature,  and  other 
thermodynamic  properties  needs  to  be  found.  This  will  allow  a 
comparison  of  experimental  data  to  theoretical  predictions. 

Experimental  data  of  differential  volume  as  a  function  of 

14 

temperature  is  available  for  carbon  disulphide.  The  equation  of 
state  used  for  this  comparison  is  Van  der  Waals'  equation  expanded 
out  to  four  terms.15  It  is  given  as 


p  +  ^  =  [1+v  +  -625(v)2  +  2869  (V)3  +  1928(v)4]- 


Temperature  as  a  function  of  volume  was  plotted  and  is  shown 
in  Figure  1.6.  The  constants  a  and  b  in  the  operating  state  were 
adjusted  to  fit  this  data.  The  adjusted  values  of  a  and  b  are  41.6 
L^-atm/mol^  and  .173  L/mol,  respectively.  The  accepted  values  of  a 
and  b  for  CS2  gas  are  given  as  11.62  L^atm/mol^  and  .07685  L/mol, 
respectively. 16  Our  values  should  be  considered  to  be  empirical 
parameters  only. 
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Figure  1.6 


References 


1.  P.  J.  Westervelt  and  R.  S.  Larson,  "Laser-excited  Broadside 
Array,"  J.  Acoust.  Soc.  Am.  £4,  121  (1973). 

2.  D.  A.  Hutchins  and  A.  C.  Tam,  "Pulsed  Photoacoustic 
Materials  Characterization,"  IEEE  Transactions  on  Ultrasonics, 
Ferroelectronics  and  Frequency  Control,  vol.  UFFC-33  No. 5  Sept. 
1986. 

3.  A.  C.  Tam,  "Applications  of  Photoacoustic  Sensing  Techniques," 
Rev  of  Modern  Physics,  Vol.  58,  No.2,  381,  April  1986. 

4.  Y.  H.  Berthelot  and  I.  J.  Busch- Vishniac,  "Laser-induced 
Thermoacoustic  Radiation,"  J.  Acoust.  Soc.  Am.  2ft,  2074  (1984). 

5.  B.  Sullivan  and  A.  C.  Tam,  "Profile  of  Laser-produced  Acoustic 
Pulse  in  a  Liquid,"  J.  Acoust.  Soc.  Am.Zft,  437  (1984). 

6.  H.  M.  Lai  and  K.  Young,  "Theory  of  the  Pulsed  Optoacoustic 
Technique,"  J.  Acoust.  Soc.  Am.22,  2000  (1982). 

7.  J.  H.  Andreae,  E.  L.  Heasel  and  J.  Lamis,  "Ultrasonic  Relaxation 
and  the  Vibrational  Specific  Heat  of  Carbon  Disulfide,"  Proc. 
Phys.  Soc.  B69,  625  (1956). 

8.  Charles  H.  Thompson,  Stanley  A.  Cheyne,  H.  E.  Bass  and  Richard 
Raspet,  "Optoacoustic  Observation  of  Internal  Relaxation  in 
Liquid  CS2,"  submitted  to  J.  Acoust.  Soc.  Am.,  8/88. 

9.  Y.  Higashigaki,  S.  L.  Whittenburg,  and  C.  H.  Wang,  "Light 
Scattering  Studies  of  Orientational  Fluctuations  of  CS2,"  J.  Chem. 

Phys.  fi2(7),  3297(1978). 

10.  G.  S.  Bushanam,  F.  S.  Barnes,  "Laser-generated  Thermoelastic 
Shock  Wave  in  Liquids,"  J.  Appl.  Phys.  46,  2074  (1975). 

11.  Handbook  of  Chemistry  and  Physics.  61st  Edition,  CRC  Press, 
p.C-231. 


20 


12.  Handbook  of  Chemistry  and  Physics. 61st  Edition,  CRC  Press, 
p.D-176. 

13.  Handbook  of  Chemistry  and  Phvsics.61st  Edition,  CRC  Press, 
p.C-682 

14.  American  Institute  of  Physics  Handbook,  edited  by  D.  E.  Gray, 
(McGraw-Hill,  1972),  p2-178 

15.  Hirshfelder,  Curtiss,  and  Bird,  Molecular  Theory  of  Gases  and 
Liquids.  (John  Wiley  &  Sons,  Inc.,  New  York,  London),  p.282. 

16.  Handbook  of  Chemistry  and  Physics.  61st  Edition,  CRC  Press,  p. 
D-194. 


21 


2.0 


During  the  previous  one  year  period,  this  lower  prior  task  has 
not  been  actively  pursued  due  to  difficulties  in  obtaining  polymer 
samples  of  interest. 
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3.0  Propagation  of  Sound  in  Vibrationally  Excited  Gases 

(SACER) 


3.1  Summary  of  Past  Work 

This  project  was  begun  several  years  ago  as  a  search  for  SACER 
(Sound  Amplification  from  Controlled  Excitation  Reactions).  The 
effect  has  been  observed  in  N2  with  added  amounts  of  He,  H2,  H20 
and  CH4  and  the  measured  amplification  found  to  exceed 
theoretically  predicted  values.  In  the  process,  the  experimental 
procedure  has  enabled  the  measurement  of  some  here-to-fore 
difficult-to-observe  relaxation  times  and  their  temperature 
dependence.  These  results  were  reported  in  a  paper  entitled 
"Propagation  of  Sound  in  Vibrationally  Excited  N2/He,  N2/H2, 
N2/CH4,  and  N2/H20  Mixtures"  published  in  JASA  &1,  2186  (1988). 

The  experimental  method  has  been  described  in  earlier  reports. 
It  involves  exciting  the  vibrational  states  in  a  gas  with  an  electrical 
discharge  and  observing  the  change  in  velocity  and  absorption  of 
sound  in  the  gas  over  a  period  of  a  few  tens  of  milliseconds  following 
the  discharge.  In  N2  and  a  few  other  selected  gases  most  of  the 
electrical  discharge  energy  winds  up  in  molecular  vibrational  energy 
in  a  few  microseconds.  By  measuring  the  energy  the  discharge 
dumps  into  the  gas,  we  have  estimated  the  vibrational  temperature 
immediately  following  the  discharge.  Vibrational  temperatures  1000 
to  2000  K  above  translational  temperatures  are  obtained  in  this  way. 
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Translational  temperatures  are  determined  from  the  sound 
velocity  which,  in  turn,  is  determined  from  the  resonant  frequency  of 
the  discharge  tube  which  changes  with  time  following  the  discharge. 
The  sound  velocity  in  this  way  becomes  a  thermometer  able  to  sense 
the  gas  temperature  changes  on  a  millisecond  time  scale,  and  the 
discharge  tube  a  kind  of  spectraphone  that  measures  the  rate  of 
conversion  of  vibrational  energy  into  translational  energy. 

3.2  Work  During  the  Past  Year 

We  have  spent  the  past  year  working  primarily  with  carbon 
monoxide.  This  gas  was  chosen  because  calculations  from  electron 
collision  cross  sections  showed  that  it  could  also  be  excited 
vibrationally  by  an  electrical  discharge.1  It  was  expected  to  behave 
in  a  way  similar  to  N2  but  to  have  the  additional  property  of  being 
infrared  active.  We,  therefore,  expected  to  be  able  to  observe  the 
decay  of  the  vibrational  temperature  directly  from  the  emitted 
infrared.  (In  N2  this  decay  was  only  inferred  from  the  resulting  rise 
and  fall  of  the  translational  temperature).  Tim  Ruppel  is  working  on 
this  problem  for  his  Master's  Thesis.  Efforts  to  date  can  be 
summarized  under  three  headings: 

1 )  Bibliography  of  previous  work  with  CO, 

2)  Efforts  at  vibrationally  exciting  CO  with  an  electrical 
discharge, 

3)  Measurement  of  vibrational  temperatures  in  N2  using  CO  as 


the  monitoring  gas. 

These  three  areas  are  discussed  separately  below. 
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3.2b.  Efforts  at  Vibrationally  Exciting  CO  in  an  Electrical  Discharge 


So  far  we  have  been  successful  in  getting  only  small  amounts  of 
energy  into  CO  and  CO  light  gas  mixtures  with  our  electrical 
discharge.  This  is  evidently  because  of  the  larger  conductivity  of  the 
CO  gas.  We  are  continuing  an  effort  to  find  a  discharge  mechanism 
that  give  us  greater  vibrational  excitation  of  CO. 

3.2c.  Measurement  of  the  Vibrational  Temperature  of  N2  Using 
CO  as  a  Monitoring  Gas 

Perhaps  the  most  significant  accomplishment  of  the  past  year 
has  been  the  successful  measurement  of  vibrational  temperatures  in 
N2  using  CO  as  an  indicator  gas.  The  vibrational  energy  levels  in  CO 
and  N2  are  sufficiently  close  so  that,  when  a  trace  amount  of  CO  is 
mixed  with  N2,  the  CO  vibrational  energy  will  maintain  thermal 
equilibrium  with  the  N2  vibration  and  the  two  vibrational 
temperatures  decay  together.2  The  infrared  activity  of  CO  can 
therefore  be  used  to  follow  the  vibrational  energy  of  N2  even 
through  the  non-polar  N2  molecule  does  not  radiate. 

The  vibrational  temperatures  have  been  measured  using  a  cold 
gas  cell  method.3  For  this  purpose,  the  electrical  discharge  tube  has 
been  terminated  with  a  cylindrical  cell  7.31  cm  long  with  CaF  end 
windows.  (See  Fig.3.1)  When  CO  gas  at  atmospheric  pressure  is 
introduced  into  this  cell  at  room  temperature,  the  part  of  the 
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infrared  corresponding  to  the  1  to  0  vibrational  transition  will  be 
absorbed.  The  ratio  of  the  infrared  intensity  coming  through  the 


empty  cell  to  that  coming  through  the  full  cell  should  be  given  by  the 


equation 


I  (empty  )  -I  (full)  \-2 

I  (  empty )  1  ' 


Where  e  =  exp  [-h-^  -1 
kTvJ 


This  equation  therefore,  can  be  solved  to  give  the  vibrational 
temperature. 

Figure  3.2  shows  results  of  measurements  on  a  mixture  of  N2 
and  CO.  Curve  1  is  the  translational  temperature  as  determined  from 
the  sound  velocity  measurements.  The  translational  temperature 
changes  little  with  time  due  to  the  slow  v-t  rate  in  pure  N2.  Curve  2 
is  the  vibrational  temperature  of  N2  calculated  from  the  translational 
temperature  assuming  a  relaxation  model  and  known  values  for  the 
thermal  conduction  and  vibrational  relaxation  times  as  explained  in 
ref.4.  Curve  3  is  the  vibrational  temperature  of  CO  as  determined 
from  the  decay  of  the  CO  infrared  radiation  using  the  "cold-gas-cell" 
method  described  above.  The  difference  between  curves  2  and  3 
caused  us  no  little  concern  until  we  remembered  that  even  through 
the  vibrational  modes  in  N2  and  CO  are  in  equilibrium,  the  principle 
of  detailed  balancing  predicts  that  their  vibrational  temperatures 
will  not  be  equal  if  the  vibrational  temperatures  differ  from  the 
translational  temperature.  The  relationship  between  the  vibrational 
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temperature  of  N2  (TN2),  the  vibrational  temperature  of  CO  (Tco),  and 
the  translational  temperature  (T)  is  given  by3 


Vco  -  VN2  _  Vco  _  VN2 
T  Tco  Tn2 

where  vco  and  VN2  are  the  fundamental  vibrational  frequencies  of 
the  CO  and  N2  molecules  using  this  equation,  curve  4  in  Fig.3.2  is 
calculated  from  curves  2  and  1.  This  agreement  between  the 
measured  relaxation  of  the  vibration  temperature  (curve  4)  and  that 
inferred  from  the  energy  deposited  in  the  gas  by  the  electrical 
discharge  (curve  2)  is  very  gratifying  and  gives  us  confidence  in  our 
method  of  determining  relaxation  times  and  their  temperature 
dependence. 


3.3  Work  Planned  for  the  Coming  Year 

The  discovery  that  we  can  use  CO  as  an  "indicator  gas"  to 
measure  the  decay  of  vibrational  temperatures  in  our  SACER- 
SPECTRAPHONE  opens  up  some  interesting  possibilities.  One  that  we 
want  to  pursue  during  the  coming  year  is  the  measurement  of  the 
vibrational  relaxation  of  H2.  The  very  large  spacing  between  the 
vibrational  energy  levels  in  H2  means  that  H2  gas  must  be  at  several 
thousand  degrees  before  the  gas  contains  appreciable  vibrational 
energy.  For  that  reason  no  one  has  ever  measured  vibrational 
relaxation  rates  in  H2.  Electron  collision  cross  section  values  indicate 
that  H2  like  N2  and  CO  should  be  vibrationally  excited  by  an  electric 
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discharge.  In  addition,  the  second  vibrational  level  in  CO  is  matched 
almost  exactly  with  the  first  level  in  H2.  Therefore  v-v  exchange 
between  CO  ahd  H2  should  be  fast  and  we  should  be  able  to  use  CO  as 
an  "indicator  gas"  to  follow  the  decay  of  vibrational  energy  in  H2  as 
was  discussed  above  for  N2.  If  this  is  possible  the  variation  of  the 
relaxation  time  in  H2/CO  mixtures  with  CO  concentration  should  give 
some  valuable  information  not  previously  available  about  v-t  and  v- 
v  rates  in  this  and  other  H2  gas  mixtures. 
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